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ABSTRACT 

Using data obtained by the EUV Imaging Spectrometer on board Hinode , we have performed a survey of obvious 
and persistent (without significant damping) Doppler shift oscillations in the corona. We have found mainly two 
types of oscillations from February to April in 2007. One type is found at loop footpoint regions, with a dominant 
period around 10 minutes. They are characterized by coherent behavior of all line parameters (line intensity, Doppler 
shift, line width, and profile asymmetry), and apparent blueshift and blueward asymmetry throughout almost the 
entire duration. Such oscillations are likely to be signatures of quasi-periodic upflows (small-scale jets, or coronal 
counterpart of type-II spicules), which may play an important role in the supply of mass and energy to the hot 
corona. The other type of oscillation is usually associated with the upper part of loops. They are most clearly seen 
in the Doppler shift of coronal lines with formation temperatures between one and two million degrees. The global 
wavelets of these oscillations usually peak sharply around a period in the range of three to six minutes. No obvious 
profile asymmetry is found and the variation of the line width is typically very small. The intensity variation is often 
less than 2%. These oscillations are more likely to be signatures of kink/ Alfven waves rather than flows. In a few 
cases, there seems to be a n / 2 phase shift between the intensity and Doppler shift oscillations, which may suggest 
the presence of slow-mode standing waves according to wave theories. However, we demonstrate that such a phase 
shift could also be produced by loops moving into and out of a spatial pixel as a result of Alfvenic oscillations. 
In this scenario, the intensity oscillations associated with Alfvenic waves are caused by loop displacement rather 
than density change. These coronal waves may be used to investigate properties of the coronal plasma and 
magnetic field. 
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1. INTRODUCTION 

In the past 15 years, there was a rapid growth of literature 
reporting quasi-periodic oscillations in the upper solar atmo- 
sphere (for latest reviews, see, e.g., Banerjee et al. 2007; De 
Moortel & Nakariakov 2012). These oscillations are often ex- 
plained as signatures of various modes of MHD waves. For 
example, Jess et al. (2009) reported an oscillation associated 
with a large magnetic bright point (BP) and interpreted it as 
a torsional Alfven wave. Intensity oscillations of a flare loop 
observed in radio wavelengths were interpreted as fast sausage 
mode and applied to constrain the density contrast between the 
loop and the background (Nakariakov et al. 2003; Aschwanden 
et al. 2004). Evidence of fast sausage modes was also recently 
found in magnetic pores (Morton et al. 2011). 

Lateral transverse oscillations are often interpreted as fast- 
mode kink waves or Alfven waves. The ubiquitous presence of 
transverse waves has been observed in chromospheric spicules 
and mottles (De Pontieu et al. 2007; Singh & Dwivedi 2007; 
Zaqarashvili et al. 2007; Kukhianidze et al. 2006; Kim et al. 
2008; He et al. 2009b; Okamoto & De Pontieu 2011; Kuridze 
et al. 2012), coronal loops (Aschwanden et al. 1999; Nakariakov 
et al. 1999; Nakariakov & Ofman 2001; Verwichte et al. 2004; 
Tomczyk et al. 2007; Tomczyk & McIntosh 2009; McIntosh 
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et al. 2011; Singh et al. 2011), and plume-like structures 
in different regions of the Sun (McIntosh et al. 2011; Tian 
et al. 2011c). Case studies suggest that some transverse waves 
are likely associated with impulsive phenomena such as flare 
activities, filament eruptions, or coronal mass ejections (CME; 
Aschwanden et al. 1999, 2002; Schrijver et al. 2002; Wang 
& Solanki 2004; O'Shea et al. 2007b; Chen et al. 2008; 
Ofman & Wang 2008; He et al. 2009a; Verwichte et al. 2010a; 
Aschwanden & Schrijver 2011; Pietarila et al. 2011; Liu et al. 
2011, 2012b; Shen & Liu 2012; White & Verwichte 2012; White 
et al. 2012; Wang et al. 2012b). Large-scale transverse waves 
with clear signatures of damping have also been reported in 
coronal streamers and found to be driven by CME impingement 
(Chen et al. 2010, 201 1 ; Feng et al. 2011). Modeling efforts have 
also been made to understand the generation and propagation 
of these transverse waves (e.g., Matsumoto & Shibata 2010; 
Murawski & Musielak 2010a; Soler et al. 2011; Selwa & Ofman 
2010; Selwa et al. 2011a, 2011b) and we refer to a recent 
review of Ofman (2009). Observations of these Alfvenic waves 
(Goossens et al. 2009), either persistent or impulsively excited, 
have been widely used for the derivation of coronal Alfven speed 
and magnetic field. For reviews of the observations of these 
transverse waves and their applications in coronal seismology, 
we refer to Nakariakov & Verwichte (2005), Erdelyi (2008), and 
Zaqarashvili et al. (2009). In spite of the implications of such 
transverse waves and their role in the MHD seismology of the 
solar corona, the detection of such waves may also be important 
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for diagnosing the stellar coronae (e.g., Mitra-Kraev et al. 2005; 
Pandey & Srivastava 2009). 

Longitudinal oscillations have also been intensively studied. 
Using mainly broadband imaging observations, quasi-periodic 
upward propagating disturbances with a period of 3-30 minutes 
were frequently observed in polar plumes (e.g., Ofman et al. 
1997, 1999, 2000; DeForest & Gurman 1998; Krishna Prasad 
et al. 2011) and coronal loops (e.g., Berghmans & Clette 
1999; De Moortel et al. 2000, 2002; Robbrecht et al. 2001; 
Marsh et al. 2003, 2009; King et al. 2003; De Pontieu et al. 
2003; Ugarte-Urra et al. 2004a, 2004b; De Pontieu et al. 2005; 
McEwan & De Moortel 2006; Tian & Xia 2008; Srivastava 
et al. 2008; Srivastava & Dwivedi 2010; Stenborg et al. 2011; 
Yuan & Nakariakov 2012; Krishna Prasad et al. 2012b). A 
propagating disturbance with a shorter period (50 s) has been 
recently reported by Liu et al. (2012a) in a dark cavity region. 
Longitudinal disturbances have also been investigated through 
spectroscopic observations (e.g., Banerjee et al. 2000, 2009; 
O'Shea et al. 2007a; Gupta et al. 2009, 2010; Wang et al. 
2012a; Krishna Prasad et al. 2012a), which often reveal a 
correlation between changes of the line intensity and Doppler 
shift (Brynildsen et al. 1999; Wikstpl et al. 2000; Wang et al. 
2009a, 2009b; Kitagawa et al. 2010; Mariska & Muglach 2010; 
Gupta et al. 2012). These long-lasting propagating disturbances 
often have a speed of 50-200 km s -1 and are usually interpreted 
as slow-mode magnetoacoustic waves propagating into the 
corona along the magnetic field lines. De Pontieu et al. (2004) 
proposed that some of these longitudinal oscillations might be 
caused by the leakage of p-modes into higher layers of the solar 
atmosphere along inclined magnetic field lines. Longitudinal 
oscillations with clear signatures of damping have also been 
detected (Wang et al. 2002, 2003b, 2003a, 2005; Mariska et al. 
2008). In some cases, there is a 7t/2 phase shift between the 
oscillations of the intensity and Doppler shift (Wang et al. 2003b, 
2003a; Mariska et al. 2008), suggesting the possible presence 
of standing slow waves. There are suggestions that these 
longitudinal oscillations can be used to infer thermal conduction 
(Ofman & Wang 2002), adiabatic index (Van Doorsselaere 
et al. 2011), and magnetic field (Wang et al. 2007) in coronal 
loops. For reviews of the observations and theories of these 
long-lasting longitudinal propagating disturbances and quickly 
damped longitudinal oscillations, we refer to De Moortel (2009), 
Taroyan & Erdelyi (2009), Banerjee et al. (2011), and Wang 
( 2011 ). 

However, De Pontieu & McIntosh (2010) and Tian et al. 
(2011a) argued that quasi-periodic signals are not necessarily 
waves. They found that in loop footpoint regions, all of the 
four line parameters (line intensity, Doppler shift, line width, 
and profile asymmetry) show coherent behavior, which can 
be easily explained by the scenario of recurring upflows (or 
quasi-periodically enhanced and weakened upflows). In addi- 
tion, Tian et al. (2011a) found a net blueshift and blueward 
asymmetry of the line profiles at the boundary of an active re- 
gion (AR) throughout the entire observation duration, favoring 
the interpretation of a recurring high-speed (50-150 km s -1 ) 
outflow superimposed on the nearly static coronal background. 
The presence of these upflows has been confirmed through de- 
tailed analysis of emission line profiles (Hara et al. 2008; De 
Pontieu et al. 2009; McIntosh & De Pontieu 2009a, 2009b; Peter 
2010; Bryans et al. 2010; Martinez- Sykora et al. 2011; Dolla & 
Zhukov 2011; Ugarte-Urra & Warren 2011; Tian et al. 2011a, 
201 lb, 2012; Doschek 2012; Rouppe van der Voort et al. 2009). 
Tian et al. (2012) also reported similar high-speed upflows in 
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CME-induced dimming regions. Based on the similarities of 
the upward propagating disturbances in plumes and AR bound- 
aries, McIntosh et al. (2010b) and Tian et al. (2011c) suggested 
that the propagating disturbances in plumes might also be dom- 
inated by flows, although an accurate investigation of the line 
profile asymmetry in coronal holes is impossible with current 
instruments. If this scenario is correct, then these high-speed 
episodic upflows/jets might be an efficient means to provide 
heated mass into the corona and solar wind (Sakao et al. 2007; 
De Pontieu et al. 2009; McIntosh & De Pontieu 2009b; Hansteen 
et al. 2010; De Pontieu et al. 2011; Tian et al. 201 lc). Recurrent 
plasma ejections and X-ray jets have also been identified from 
recent imaging observations (Morton et al. 2012; Zhang et al. 
2012; Su et al. 2012; Xia et al. 2005). Several numerical stud- 
ies have been performed to investigate the formation of these 
plasma ejecta (e.g., Murawski & Zaqarashvili 2010b; Murawski 
et al. 2011; Srivastava & Murawski 2011; Martinez- Sykora et al. 
2011; McLaughlin et al. 2012). 

The discovery of these rapid quasi-periodic upflows (or 
jets, plasma ejections) from spectroscopic observations chal- 
lenges the universal wave interpretation for coronal oscillations, 
specifically for some longitudinal propagating disturbances. 
There is no doubt that both waves and flows exist in the up- 
per solar atmosphere, but it is difficult to distinguish between 
them without analyzing all line parameters including intensity, 
Doppler shift, line width, and profile asymmetry (De Pontieu & 
McIntosh 2010; Tian et al. 201 la), although realistic numerical 
MHD models may help resolve this ambiguity (Ofman et al. 
2012). In this paper, we perform a survey of obvious and per- 
sistent (without significant damping) Doppler shift oscillations 
in the corona by using the data obtained by the EUV Imag- 
ing Spectrometer (EIS; Culhane et al. 2007) on board Hinode. 
Through detailed analysis of the temporal evolution of all the 
line parameters, we conclude that some oscillations are better 
explained by recurring upflows and others are more likely to be 
real waves. 

2. DATA SELECTION AND REDUCTION 

EIS made a large number of sit-and-stare (fixed slit location) 
observations from February to April in 2007. We examined the 
quick-look Dopplergrams (automatically generated level 2 data) 
of the Fe xn A 195. 12 line, which are available on the Web site of 
Hinode Science Data Center Europe, for these observations. In 
Table 1 , we list the starting time of all observations in which we 
found obvious (easily identified by eye) and persistent (without 
significant damping) Doppler shift oscillations. The locations 
of oscillations on the slits are also listed in Table 1 . The typical 
cadence of these observations is 30 s. 

After standard correction and calibration of the EIS data, a 
running average over three pixels along the slit was applied to the 
spectra to improve the signal-to-noise ratio. The spatial offset in 
the solar T-direction between the two CCDs has been corrected 
using the standard EIS routine eis_ccd_offset.pro. This routine 
measures the offset by co-aligning images from Fe viii A 185. 21 
and Si vn A275.35. We mainly selected three strong emis- 
sion lines in each observation for our study: Fexn A 195. 12, 
Fexm A202.04, and Fexiv A264.78. The latter two lines are 
believed to be clean (Young et al. 2007), although an uniden- 
tified line might be present in the red wing of Fe xiii A202.04 
(McIntosh et al. 2010a; Tian etal. 2012). The Fe xn A 195. 12 line 
is known to be blended with the Fe xii A 1 95 . 1 8 line that is usually 
a few percent of the 195. 12 A intensity (Young et al. 2007, 2009; 
Brown et al. 2008). We used this line since it is the strongest 
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Table 1 

Persistent Doppler Shift Oscillation Events Observed by EIS from February to April in 2007 


id 

Date 

Time 

Y 

Ri 

Rv 

Rw 

Ra 

Pi 

Piun 

Pv 

Pvun 

Ai 

Av 

N 

Nerr 

Type, Location 

1 

20070201 

013212 

351 

1.07 

0.44 

0.30 

0.46 

10.7 

2.44 

9.82 

2.05 

3.89 

1.65 

9.23 

0.04 

I, loop leg 

2 

20070202 

004912 

245 

0.89 

0.49 

0.51 

0.65 

7.57 

3.19 

10.7 

4.79 

2.36 

1.45 

9.25 

0.04 

I, loop leg 

3 

20070203 

005642 

193 

1.14 

0.59 

0.56 

0.63 

8.77 

3.00 

7.57 

2.44 

2.71 

1.60 

9.22 

0.04 

I, loop leg 

4 

20070220 

175013 

296 

4.57 

2.79 

2.33 

2.60 

11.6 

2.04 

12.7 

2.45 

14.0 

8.47 

9.65 

0.13 

I, loop leg 

5 

20070221 

021812 

310 

3.45 

1.62 

1.33 

3.12 

10.7 

2.90 

9.82 

2.05 

11.8 

5.19 

9.32 

0.17 

I, loop leg 

6 

20070326 

150012 

168 

0.87 

0.48 

0.48 

0.51 

8.26 

1.88 

9.82 

2.05 

2.53 

1.45 

9.21 

0.05 

I, loop leg 

7 

20070327 

143412 

220 

1.24 

0.56 

0.41 

0.49 

13.8 

2.42 

9.82 

4.89 

4.11 

1.40 

9.65 

0.10 

I, loop leg 

8 

20070328 

031531 

264 

1.84 

0.89 

0.77 

0.81 

13.8 

3.17 

11.7 

2.99 

5.50 

2.63 

9.03 

0.05 

I, loop leg 

9 

20070328 

183726 

303 

2.08 

0.56 

0.77 

0.76 

9.01 

2.64 

10.7 

2.70 

6.25 

1.67 

9.52 

0.07 

I, loop leg 

10 

20070216 

163920 

345 

1.93 

0.70 

0.46 

0.71 

5.35 

1.29 

4.13 

1.46 

4.15 

1.43 

9.67 

0.07 

II, QS BP 

11 

20070221 

194813 

394 

1.22 

0.54 

0.34 

0.56 

10.7 

1.34 

3.47 

0.43 

4.73 

1.57 

8.83 

0.05 

II, loop top 

12 

20070223 

125943 

384 

1.09 

0.46 

0.22 

0.45 

6.37 

1.46 

4.91 

2.19 

2.09 

0.95 

9.27 

0.07 

II, loop top 

13 

20070223 

185342 

380 

1.06 

0.56 

0.35 

0.55 

6.14 

0.87 

3.47 

0.61 

2.38 

1.15 

9.25 

0.08 

II, loop top 

14 

20070224 

005742 

381 

1.27 

0.54 

0.34 

0.57 

5.84 

1.68 

3.18 

0.67 

2.48 

1.08 

9.31 

0.09 

II, loop top 

15 

20070224 

130643 

370 

0.88 

0.63 

0.25 

0.31 

4.50 

1.03 

3.47 

0.30 

1.67 

1.76 

9.32 

0.09 

II, loop top 

16 

20070225 

112113 

185 

0.91 

0.66 

0.29 

0.46 

7.57 

3.06 

4.91 

1.29 

1.95 

1.58 

9.10 

0.05 

II, unclear 

17 

20070326 

150012 

29 

0.79 

0.68 

0.27 

0.38 

8.26 

2.92 

5.35 

0.93 

2.21 

2.13 

9.00 

0.04 

II, loop top 

18 

20070327 

160926 

247 

0.80 

0.74 

0.28 

0.38 

6.94 

2.05 

4.91 

1.12 

2.04 

1.70 

9.16 

0.03 

II, loop top 

19 

20070328 

031531 

12 

0.84 

0.72 

0.27 

0.46 

8.26 

3.44 

7.57 

2.59 

2.03 

1.76 

8.94 

0.04 

II, loop top 

20 

20070328 

183726 

40 

0.64 

0.67 

0.26 

0.33 

12.7 

1.01 

4.91 

0.67 

1.78 

1.77 

9.14 

0.03 

II, loop top 

21 

20070419 

191102 

189 

3.17 

0.97 

0.44 

0.65 

6.94 

1.02 

5.35 

1.29 

8.39 

2.73 

9.08 

0.03 

II, QS BP 


Notes. The following information is listed for each event: observation date (yyyymmdd), starting time (hhmmss), lowest y-pixel of the selected 1 1-pixel region 
on the slit (F), root-mean-square values of the detrended intensity (Ri, %), Doppler shift (Rv, km s -1 ), line width (Rw, km s -1 ) and RB asymmetry (Ra, %), 
intensity oscillation period and uncertainty (Pi and Piun, minutes), Doppler shift oscillation period and uncertainty (Pv and Pvun, minutes), amplitudes of 
intensity (Ai, %) and Doppler shift (Av, km s -1 ) oscillations, electron density, and uncertainty (N and Nerr, log cm -3 ), and oscillation type and location. 


EIS coronal line in non-flare conditions. The Fe xn X 195. 12 line 
was also involved in density diagnostics and the co-alignment 
between Transition Region and Coronal Explorer (TRACE) and 
EIS. More lines were selected for a comprehensive analysis of 
the 2007 March 28 observation. 

As a common practice, we applied a single Gaussian fit to 
each EIS spectrum to derive the line intensity, Doppler shift, and 
line width. Since there are no cool chromospheric lines in the 
EIS spectral range, we simply assumed a zero average Doppler 
shift of each line in each observation. In some observations, 
there was significant jitter in the solar y-direction. Similar 
to Kitagawa et al. (2010), we performed a cross-correlation 
between intensities of Fexn A 195. 12 at different exposures to 
evaluate and remove the jitter in the y-direction. It turned out that 
the jitters were mostly smaller than 2". Since the oscillations we 
found are associated with specific locations on the slit, we can 
safely conclude that they are not caused by possible jitters in 
the v-direction. We also calculated the average Red-blue (RB) 
asymmetry in the velocity range of 60-120 km s -1 for each 
line profile. The RB asymmetry technique was first introduced 
by De Pontieu et al. (2009) and later modified by Tian et al. 
(2011b). The RB asymmetry here refers to RB P , where the 
spectral position corresponding to the maximum intensity is 
used as the line centroid (see details in Tian et al. 201 lb). 

3. 2007 MARCH 28 OSCILLATIONS: TWO TYPES OF 
OSCILLATIONS OBSERVED ON THE SAME SLIT 

Starting from 18:37 on 2007 March 28, EIS performed sit- 
and-stare observation in an AR continuously for more than 
3 hours. Figure 1(A) shows an image of this AR in the TRACE 
195 A passband and the location of the EIS slit as determined by 
cross-correlating the EIS Fe xn A 195. 12 line intensity along the 
slit and the TRACE intensity at different v-locations. From the 
temporal evolution of the Doppler shift of Fe xiii A202.04 shown 


in Figure 1(B), we can clearly see two locations (indicated by 
the two arrows) with obvious oscillations throughout the entire 
duration. One oscillation (around y = 162") is characterized 
by some inclined elongated features of enhanced blueshift. This 
oscillation appears to be associated with the footpoints (or lower 
parts) of AR loops, as can be seen from Figure 1(A). We call it 
oscillation 1 in the following discussion. The other oscillation 
(around y = —95") is dominated by small redshift and we can 
clearly see the quasi-periodic appearance of vertical features 
with enhanced redshift. From Figure 1(A), we can see that this 
oscillation (oscillation 2) appears to be associated with the tops 
(or upper parts) of a bundle of AR loops. The slit is found 
to be almost parallel to the loop plane and cover at least 45" 
of the loop top. Note that although the edge the bright region 
around y = —70" is caught by the slit, this part does not show 
significant oscillation and it is not analyzed in the following. 

We have performed plasma diagnostics for these two regions, 
revealing clear Doppler shift oscillations. As an example, 
Figures 1(C) and (D) show the electron densities and differential 
emission measure (DEM) curves at the positions of y = 162" 
and y = —95". The densities were derived from the intensity 
ratio of the line pair Fexn AA196. 64/195. 12. The theoretical 
relationship between the line ratio and density, as extracted 
from the CHIANTI database (Dere et al. 1997; Landi et al. 
2006), is shown as the black line. The electron densities at the 
two locations are log (N t / cm -3 ) = 9.54 ± 0.12 and 9.20 d= 0.04, 
respectively. More lines were used for the derivation of DEM 
curves. We averaged the profiles of each line acquired at different 
exposures and assumed an uncertainly of 10% for the intensities 
(Tripathi et al. 2008). By using the routine chianti _dem.pro (also 
used by Lee et al. 2011; Tian et al. 2012) available in SolarSoft 
(SSW) and assuming a constant pressure of 10 16 cm -3 K, 
we obtained the DEM curves at the two locations. The peak 
temperatures at the two locations are around log(T /K) = 5.85 
and 6.15, respectively. 
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( A ) TRACE 195 intensity (arbitrary unit) ( B ) Doppler shift of Fe XIII 202.04 (km/s) 



200 250 300 350 400 450 500 0 50 100 150 

X (arcsec) Minutes after 18:37 on 2007/03/28 



0.02 0.04 0.06 0.08 0.10 0.12 4.5 5.0 5.5 6.0 6.5 

Ratio of Fe XII 196.64/(195.12+195.18) log (T/K) 


Figure 1 . (A) TRACE 195 A image taken at 20:34 on 2007 March 28. The black vertical line indicates the location of the EIS slit. The arrows point to the approximate 
locations of two persistent oscillations. (B) Temporal evolution of the Doppler shift of Fexm A.202.04 along the slit. The short horizontal bars in (A) and dashed lines 
in (B) are used to indicate the ranges of illustration in Figures 2-5. (C) Electron densities at the positions of y = 162" and y = —95", respectively. Error bars indicate 
the standard deviations of the density time series. (D) DEM curves at the positions of y = 162" and y = —95", respectively. 

(A color version of this figure is available in the online journal.) 


3.1. Oscillation 1: Correlated Changes of All Line Parameters 

These oscillations are also found in several other emission 
lines with different formation temperatures. From the left panels 
of Figure 2, we can see that oscillation 1 is clearly present in all 
lines with a formation temperature of 0. 6-2.0 MK. The inclined 
elongated features of enhanced blueshift are most clearly seen 
in the strongest line Fexn A. 195. 12 and can also be identified 
in other weaker lines. For all coronal lines, we see blueshift 
throughout almost the entire observation duration. The transition 
region line Sivu A275.35 reveals a mixture of both blueshift 
and redshift, which might be related to the enhancement of 
downflows at lower temperatures (Tian et al. 2011b; McIntosh 
et al. 2012). We averaged all the line parameters (intensity, 
Doppler shift, line width, and RB asymmetry) over the 1 1 pixels 
between the two dashed lines and produced a time series of each 
line parameter. By subtracting a 10 minute running average from 
each time series (e.g., Wang et al. 2009a), we obtained a time 
series of each detrended line parameter, which is shown in the 
right panels of Figure 2. We found that the intensity variation 
generally follows the variation of the line width. It is also 
clear that the Doppler shift and RB asymmetry show coherent 
behaviors. By taking the inverted values of Doppler shift and RB 


asymmetry (now positive values mean blueshifts or blueward 
asymmetries), we found that all of the four parameters show 
coherent variations. Note that the less-than-ideal correlations at 
some instances are probably partly caused by the poor spectral 
resolution and photon noise of the EIS instrument (De Pontieu 
& McIntosh 2010; Tian et al. 2011a). As demonstrated by De 
Pontieu & McIntosh (2010), the photon noise has progressively 
worse effects as one goes to higher moments, i.e., the asymmetry 
suffers much more from noise than the intensity. The foreground 
and background emission in the line of sight (LOS) direction 
might be another reason for the reduced correlations. 

The detailed temporal evolutions of all of the four line pa- 
rameters of the strong Fexm A202.04 line are presented in 
Figures 3(A)-(D). The loop footpoint region is clearly char- 
acterized by quasi-periodic enhancement of the line inten- 
sity, blueshift, and line width, similar to what De Pontieu & 
McIntosh (2010) and Tian et al. (2011a) found. We can also see 
weak signatures of blueward asymmetry in the region between 
the two dashed lines. The Fe xiv A 264. 7 8 line is clean and from 
panel (F) we can clearly see obvious blueward asymmetry in 
the loop footpoint region. A comparison between panels (D) 
and (F) suggests that the Fexm A202.04 line is possibly af- 
fected by an unidentified blend (and/or slight gradient in the 
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Figure 2. Left: temporal evolution of the Doppler shift of several emission lines in the range of y = 130" ~ 210" in Figure 1(B). Right: the black, red, green, and violet 
curves represent, respectively, the detrended line intensity, Doppler shift (inverted), line width, and RB asymmetry (inverted) averaged over the region between the 
two dashed lines shown in the left panels. The intensity has been normalized to the local background and is shown as the percentage divided by 2. The RB asymmetry 
has been normalized to the peak intensity of the line profile and is shown in the unit of percentage. The units of Doppler shift and line width are km s -1 . For the 
purpose of illustration, the intensity and RB asymmetry are offset by 3 and —4, respectively, on the y-axis. 

(A color version of this figure is available in the online journal.) 


background emission) at its red wing (De Pontieu & McIntosh 
2010; McIntosh et al. 2010a; Tian et al. 2012). This effect tends 
to shift the RB asymmetry to the redward and thus reduce the 
magnitude of blueward asymmetry in the loop footpoint region. 
However, we still used the Fe xiii A 202. 04 line for our detailed 
analysis since it is much stronger than Fe xiv A264.78. The sim- 
ilar behavior of different coronal lines as shown in Figure 2 
suggests that the effect of the possible blend or background 
issue is not significant in our study. Panel (L) presents a typi- 
cal line profile of Fexiv A264.78 (nine profiles averaged) and 
we can clearly see an enhancement of the blue wing. Such an 
asymmetric line profile suggests at least two emission compo- 
nents (Hara et al. 2008; De Pontieu et al. 2009; De Pontieu & 
McIntosh 2010; McIntosh & De Pontieu 2009a, 2009b; Peter 
2010; Bryans et al. 2010; Martinez- Sy kora et al. 2011; Dolla & 
Zhukov 2011; Ugarte-Urra & Warren 2011; Tian et al. 2011a, 
2011b; Doschek 2012) and a study of center to limb variation 
clearly suggests that the enhancement in the blue wing is not 
caused by blends or noise (Tian et al. 2012). The correlated 
changes of different line parameters can be seen from the de- 
trended time series in both panels (E) and (G). Panels (J) and 


(K) show the cross-correlations between the Doppler shift and 
intensity /line width/RB asymmetry as a function of time lag. 
The nearly zero time lag for each parameter pair confirms the 
coherent behaviors of all line parameters. Similar to Wang et al. 
(2009a), we measure the oscillation period as the value corre- 
sponding to the peak global wavelet power and the uncertainty 
as the half-width at half-maximum. The oscillation periods of 
the intensity and Doppler shift are 9.01 zb 2.64 minutes and 
10.7 ± 2.70 minutes, respectively. The Doppler shift period of 
10.7 ±2.70 minutes is consistent with the typical time gap be- 
tween two inclined stripes in the left panels of Figure 2. The 
large uncertainties of periods suggest that the time gap between 
two events (i.e., upflows, see below) varies significantly with 
time. The oscillation amplitude, defined by the square root of 
the peak global wavelet power (Wang et al. 2009a), is 6.25% for 
the intensity and 1.67 km s -1 for the Doppler shift. Note that 
the data gaps separate the entire time series into three segments 
and that we selected the longest segment for wavelet analysis. 
We have also performed the wavelet analysis for the other two 
shorter segments and the obtained periods are close to that of 
the central segment. 
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Figure 3. (A-D) Temporal evolution of the line parameters (intensity: detrended and others: original) of Fexm 7202.04 in the range of y = 130" ~ 210" in 
Figure 1(B). (E) Temporal evolution of the detrended line intensity (i), Doppler shift (v), line width (w), and RB asymmetry, (a) averaged over the region between the 
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(A color version of this figure is available in the online journal.) 


6 


The Astrophysical Journal, 759:144 (17pp), 2012 November 10 


Tian et al. 


Doppler shift (km/s) 


-20 


-10 


10 


20 


^ -60 pFe XII 195.12 

0 
0 

0 -100E 


-80 r 


>■ -120 


-80 r 


^ -60 pFe XIII 202.04 

0 
0 

§ -loo : 


> -120 


o 
0 
0 

js -100 = 

> -120 r 


-60 

-80 


o 
0 
0 

js -ioo 

>■ -120 r 




-60 ip e XIV 274.20 
-80 7 


~ 




“1 




= 


Te XV 284.16 


— ^ 


1 

1 

1 

1 

T 1 

IMj ft . i. , • 

t T 



1 


“= 


0 


50 


100 


150 


C\J 


Intensity (% X 0.5) 

Doppler shift (km/s) 

Line width (km/s) 

RB asymmetry [60-120 km/s] (%) 



i 


i 


i i i i 


o ® 7 

2! o: " 

= -2 F I v I I ' I 

^ -6 | | | | | | | | | | | | L_ 

P 6 r 


i 


i 


i 


i 


i i i i 


i 


^ i v i r ' i v y u p v iv v i v i i'i i vi ■ 

0 -4 ^/WyvaIV^i 

^ -6L_I I I I | I II I | I I I I 



0 

-2 7 V i v r v i v r v i v v u \i v y v i v i i i i v| ' 

-4 ■ 

-6|f_i i i i i i i i J i i i i L 


1 1 1 J 

i r VV v r\rA^ 

iii-= 

J I J 1= 

1 1 1 J 

i--'AT / M/'A^ 


i i v i 
I /A/V'TvWVi 

J I L 

1 1 'a J 

vw 

V^v~rv« 

J I I L 


0 


50 100 150 

Minutes after 18:37 on 2007/03/28 


Minutes after 18:37 on 2007/03/28 

—53" in Figure 1(B). The Doppler shift and RB asymmetry values are not inverted here, which is 


Figure 4. Similar to Figure 2, but for the range of y = —133" ' 
also different from Figure 2. 

(A color version of this figure is available in the online journal.) 


3.2. Oscillation 2: Oscillation Most Clearly Seen 
in Doppler Shift 

Oscillation 2 can be clearly identified in all lines with a for- 
mation temperature ranging from 1.3 to 2.0 MK. From the left 
panels of Figure 4, we can see that oscillation 2 is characterized 
by almost vertical stripes of enhanced redshift. To the best of 
our knowledge, such obvious long-lasting oscillations in warm 
(1. 3-2.0 MK) coronal lines have never been found in obser- 
vations before. Lower-temperature lines such as Sivu A275.35 
and Fex A 184.54 are weak in this loop-top region so that we 
could not study the detailed temporal evolution of their Doppler 
shift. Similar to Figure 2, we obtained time series of the de- 
trended line parameters for each line in the region between the 
two dashed lines and present them in the right panels of Figure 4. 
It is clear that the oscillation is most prominent in Doppler shift. 
The variations of the intensity, line width, and RB asymmetry 
are all typically smaller than those of oscillation 1 . 

The temporal evolution of all of the four line parameters 
of Fexm A202.04 is detailed in Figures 5(A)-(D). Clearly, 
oscillation 2 is mainly present in the Doppler shift. No clear 
discernible variations are found in the line width and RB 
asymmetry. The Fe xiv A264.78 line shows behaviors similar to 
Fe xiii A202.04, as can be seen from panels (F) and (G). The line 
profiles in this loop-top region are typically symmetric and an 
example is presented in panel (L). The variation of the detrended 
intensity is smaller than that of oscillation 1. But from panel 
(G), we see a seeming correspondence between the intensity 
and Doppler shift variations. The correspondence can also be 
seen from Figure 4 and is more clear in hotter lines such as 
Fe xiv A274.20 and Fe xv A284.16. From the analysis of cross- 
correlation, we find a maximum correlation coefficient at the 
time lag of about 0.9 minutes between the time series of intensity 
and Doppler shift for the two lines shown in Figure 5. Such a 
time lag corresponds to a phase shift of about it /2 or 1/4 period, 
given a Doppler shift oscillation period of 4.91 zb 0.67 minutes. 


However, since the correlation coefficient is smaller than 0.3 
(which can easily be caused by random, chance correlations) 
and the dominant period (12.7 minutes ± 1.01) of the intensity 
is largely different from that of the Doppler shift, we cannot 
make any solid conclusion about it. The oscillation amplitude is 
1.78% for the intensity and 1.77 km s -1 for the Doppler shift. 
Clearly, the intensity fluctuation is much smaller than that of 
oscillation 1 , while the amplitude of the Doppler shift oscillation 
is very similar to that of oscillation 1 . 

4. STATISTICAL RESULTS AND DISCUSSION 

Table 1 lists the observational information and analysis re- 
sults, including observation time, locations on the slit, root- 
mean-square values of all line parameters, oscillation ampli- 
tudes, periods and uncertainties of the intensity and Doppler 
shift, electron densities, and locations, for all oscillation events 
we identified. These parameters were derived from an aver- 
aged time series over 11 pixels along the slit in each event. 
We carefully selected these 1 1 Y-pixels in each event and made 
sure that (1) the temporal behavior at these 11 pixels is al- 
most the same and (2) that this 1 1 -pixel region covers the major 
oscillating part. Similarly, we detrended each time series by sub- 
tracting a 10 minute running average. We mainly used the strong 
Fe xiii A 202. 04 line which is available in all observations for our 
statistical studies. The weaker but clean Fe xiv A264.78 line was 
also analyzed for a confirmation of the characteristics revealed 
by Fexm A202.04. The line pair of Fexn AA196.64/195.12 or 
Fexn AA186. 88/195. 12 was used to diagnose the electron den- 
sities. Images obtained by the TRACE (Handy et al. 1999), the 
X-ray Telescope (XRT; Golub et al. 2007), and the Extreme Ul- 
traviolet Imaging Telescope (Delaboudiniere et al. 1995) were 
used for identifying the coronal structures associating with these 
oscillations. Based on these statistics, we can subdivide the os- 
cillations we found here into two types. One type (type I) are 
mainly found in the lower part of AR loops (Section 4.1) and 
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Figure 5. Similar to Figure 3, but for the range of y = —133" ~ —53" in Figure 1(B). In panels (J) and (K), only the results of cross-correlation between the Doppler 
shifts and intensities are shown. In panel (L), only the single Gaussian fit is applied to the observed line profile. 

(A color version of this figure is available in the online journal.) 


they are probably related to recurring upflows. The other type 
of oscillations (type II) are found mostly near the top of loops, 
and the observational facts seem to suggest an interpretation of 
transverse wave. We have to mention that event 6 is very close 
to a localized bright region ( TRACE 195 A passband) in the core 
of an AR, but from the XRT image, it is clear that the oscillating 
part (the localized bright region in TRACE 195 A passband) is 
the leg (or lower part) of a hot loop system. We could not find 
any context images for event 16, so its location is unclear. 


4.1. Type I Oscillations: Recurring Upflows? 

We found nine oscillation events that reveal properties similar 
to oscillation 1 discussed in the previous section. One example 
is shown in Figure 6(A). These oscillations are all found at 
the footpoint regions (lower parts) of AR loops. They are 
characterized by coherent behaviors of all line parameters (line 
intensity, Doppler shift, line width, and profile asymmetry), 
obvious blueshift and blueward asymmetry throughout almost 
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Figure 6. (A) Doppler shift and time series of all line parameters of Fe xiv A.264.78 for a type I oscillation event. The line styles, units, and illustration methods are the 
same as in Figure 3(G). (B-D) The Doppler shift and time series of all line parameters of Fe xiii A. 202. 04 for three type II oscillation events. The line styles, units, and 
illustration methods are the same as in Figure 5(E). (E and F) Wavelet (of Doppler shift, similar to Figure 5(1)) and cross-correlations analysis (similar to Figure 5(J)) 
for the event shown in panel (D). 

(A color version of this figure is available in the online journal.) 


the entire observational duration. The time lag between the 
Doppler shift and other parameters is very close to zero and 
their correlation coefficients are typically larger than 0.3 for 
all of these events. The line widths are obviously enhanced at 
the oscillation locations. These oscillations are named type I 
oscillations in the following discussions. 


For these type I oscillations, the global wavelet power is 
often distributed over a wide range of periods. The dominant 
periods (period corresponding to the peak of a global wavelet) 
range from 7.57 to 13.8 minutes for the intensity, and 7.57 
to 12.7 minutes for the Doppler shift, with an average of 
10.5 zb 2.3 minutes for the intensity and 10.3 ±1.4 minutes 
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(A) Recurring upflow 




Figure 7. Synthetic line profiles (diamonds) and time series of the line parameters in the scenarios of recurring upflows (A) and periodically enhanced and weakened 
upflows (B). In the panels of line profiles, the green/red solid line is the single/double Gaussian fit. The two dashed red lines represent the primary (background) and 
secondary (upflow) components, respectively. The red and green vertical lines indicate the line center of the primary component and single Gaussian fit, respectively. 
The single Gaussian fit parameters (i: line intensity, v: Doppler shift, and w: line width) and RB asymmetry (a) are shown in each panel. 

(A color version of this figure is available in the online journal.) 


for the Doppler shift. The oscillation amplitudes are in the 
range 2.36%-14.0% for the intensity, and 1.40-8.47 km s -1 
for the Doppler shift. The average amplitudes of the intensity 
and Doppler shift are 5.9% and 2.8 km s _1 , respectively. In two 
events starting at 17:50:13 on February 20 (previously analyzed 
by Tian et al. 2011a; Nishizuka & Hara 2011) and 02:18:12 
on February 21, we see an intensity amplitude of about 13% 
and Doppler shift amplitude of about 7 km s _1 . The root-mean- 
square values of the line width and RB asymmetry for these two 
events are about 2 km s -1 and 3%, respectively. These values 
are much larger than those of other events. The magnitudes 
of fluctuations are likely influenced by effects such as LOS 
projection and the foreground/background emission. 

The line profiles are usually enhanced in the blue wings, 
suggesting the existence of at least two emission components. 
Compared to the single Gaussian fit, the RB -guided double 
Gaussian fit (De Pontieu & McIntosh 2010; Tian et al. 2011a, 
2011b, 2012) usually does a much better job. From Figure 3(L), 
one can also conclude that the enhanced line width and signif- 
icant blueshift derived from the single Gaussian fit are actually 
caused at least partly by the superposition of a weak high-speed 
upflow component on a strong background emission component 
(Tian et al. 2011b). 

These observational facts point to the explanation of recurring 
high-speed upflows (or quasi-periodically enhanced and weak- 
ened upflows) for type I oscillations, as previously suggested 
by De Pontieu & McIntosh (2010) and Tian et al. (2011a). 


Figure 7 presents synthetic line profiles and time series of the 
line parameters in the scenarios of recurring upflows and pe- 
riodically enhanced and weakened upflows. The background 
coronal emission component and high-speed upflow component 
are represented by the two red dashed lines. Clearly, continuous 
high-speed upflows with periodic enhancement of the flow in- 
tensity can easily explain the quasi-periodic enhancement of the 
line intensity, line width, blueshift, and blueward asymmetry. 
In real observations, the flow intensity might be enhanced dif- 
ferently in different periods (clearly seen in Tian et al. 2011a), 
which can sometimes cause a deviation of the dominant period 
derived from wavelet analysis from the real period of recurrence. 

Correlated changes between intensity and Doppler shift are 
usually interpreted as propagating slow-mode magnetoacoustic 
waves (e.g., Wang et al. 2009a, 2009b; Nishizuka & Hara 
2011). However, here we can see that such a correlation also 
exists in the case of recurring upflow. It seems that the line 
width and asymmetry are the key to distinguish between slow 
waves and upflows. Verwichte et al. (2010b) argued that slow 
waves can also cause line asymmetries when the emission line is 
averaged over an oscillation period or when a quasi-static plasma 
component in the LOS is included. However, the frequency 
doubling of the line width, as shown in their Figure 3, is not 
seen in our type I oscillations. In addition, the irregular quasi- 
periodic changes (change differently in different periods) of the 
intensity and Doppler shift are also expected in the scenario 
of recurring upflows since each repetitive upflow is likely to 


10 


The Astrophysical Journal, 759:144 (17pp), 2012 November 10 

be independent of each other. So at present, we conclude that 
the type I oscillations we identified here are more likely to 
be dominated by quasi-periodically enhanced and weakened 
upflows. Heggland et al. (2009) and McLaughlin et al. (2012) 
have recently demonstrated that such quasi-periodic outflows 
can be generated by oscillatory reconnection (McLaughlin et al. 
2009). Morton et al. (2012) also proposed that the recurrent 
plasma ejections they observed are reconnection-driven and 
similar to type-II spicules. Based on an observation of a 
supersonic blob, Srivastava et al. (2012) suggested that multi- 
temperature plasma blobs could be generated by a recurrent 
three-dimensional (3D) reconnection process via the separator 
dome below the magnetic null point, between the emerging 
flux and pre-existing field lines in the lower solar atmosphere. 
Indeed, quasi-periodic flow pulses associated with a null point 
have recently been reported (Su et al. 2012). On the other hand, 
Tripathi et al. (2012) explained the hot plasma upflows observed 
in the warm loops as evidence of chromospheric evaporation 
in quasi-static coronal loops, while in the 3D simulation of 
Martinez- Sy kora et al. (2011), the chromospheric plasma is 
ejected as result of being squeezed by the magnetic tension. 

We would like to point out that slow waves may also 
exist. The disappearance of enhanced line width and blueward 
asymmetry at higher parts of loops might be caused by LOS 
effect or low signal-to-noise ratio of EIS. But it might also 
be explained by slow waves (Nishizuka & Hara 2011). Ofman 
et al. (2012) recently found that periodically driven upflows 
can excite undamped slow waves in loops. They showed that 
the properties of quasi-periodic flows are mainly related to 
driving source, while the properties of waves are determined 
by the structures of loops. It is likely that high-speed upflows 
dominate in the footpoints (lower corona) of some loops and 
slow waves dominate higher up in the loops. More recently, 
through a detailed analysis of the temperature dependence of 
the speeds of the propagating disturbances, Kiddie et al. (2012) 
found signatures of both upflows and slow waves. On the other 
hand, slow waves associated with the leakage of /7-modes may 
well exist in the chromosphere (e.g., de Wijn et al. 2009) and 
some of them might propagate further up into the corona and 
mix with the upflows. 

These high-speed upflows in lower parts of loops may play 
an important role in the supply of mass and energy to the hot 
corona and probably also the solar wind (De Pontieu et al. 2009, 
2011). There are suggestions that the blueshifts of the order of 
30 km s -1 at the boundaries of some ARs (Marsch et al. 2004, 
2008; Del Zanna 2008; Tripathi et al. 2009; Murray et al. 2010; 
Warren et al. 201 1 ; Young et al. 2012), as derived from a single 
Gaussian fit, are signatures of the nascent slow solar wind (Sakao 
et al. 2007; Harra et al. 2008; Doschek et al. 2008; He et al. 
2010; Brooks & Warren 2011; Slemzin et al. 2012; Zangrilli & 
Poletto 2012). This might be true for the upflows along open 
field lines originating from the AR boundaries, but we have to 
point out that the speeds of these outflows are around 100 km s -1 
instead of 30 km s -1 , as revealed by the RB asymmetry analysis 
and double Gaussian fit (Hara et al. 2008; De Pontieu et al. 
2009; De Pontieu & McIntosh 2010; McIntosh & De Pontieu 
2009a, 2009b; Peter 2010; Bryans et al. 2010; Martfnez-Sykora 
et al. 20 1 1 ; Dolla & Zhukov 2011; Ugarte-Urra & Warren 2011; 
Tian et al. 2011a, 2011b, 2012). In addition, some of these 
outflows are clearly associated with legs of closed field lines 
(Tian et al. 2011b; Boutry et al. 2012) and they are followed 
by cooling downflows (Ugarte-Urra & Warren 2011; Kamio 
et al. 2011; McIntosh et al. 2012) instead of propagating into 
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the interplanetary space. These outflows may be the coronal 
counterpart of type-II spicules that could supply heated mass to 
the corona (De Pontieu et al. 2011). Taking an average speed 
of 100 km s -1 and assuming a density of log (N e /cm~ 3 ) = 8, 
the mass flux density of these outflows is estimated to be about 
3.3 x 10 -10 g cm -2 s -1 if we use a temporal filling factor of 0.2. 

4.2. Type II Oscillations: Waves? 

Twelve oscillation events are found to reveal characteristics 
different from those of type I oscillations. These oscillations are 
named type II oscillations in the following discussions. These 
oscillations are found to be associated with the upper parts or 
tops of AR loops. Two of them are actually found in quiet-Sun 
BPs, which are also believed to consist of a group of miniature 
loops. The oscillation is most prominent in Doppler shift and 
sometimes the cross -correlation analysis reveals a 7r/2 phase 
shift between the intensity and Doppler shift oscillations. The 
variations of the line intensity, line width, and RB asymmetry 
are all typically smaller than those of type I oscillations. The 
average root-mean-square values of the detrended line intensity/ 
line width/RB asymmetry are 1.9%/0.83 km s -1 /l. 11% for 
type I and 1.2%/0.31 km s~ 1 /0.48% for type II oscillations. 

One type II oscillation event (event 17 in Table 1) presented 
in Figure 6(B) is characterized by almost vertical stripes of 
enhanced redshift. In addition, the slit is found to be almost 
parallel to the loop plane and covered at least 50" (coherence 
scale) of the upper segment of the loop. Such characteristics are 
very similar to the events 20 (see Figures 1 and 5) and 19 (not 
shown here), and the coherence scales are around 50". There are 
signatures of correspondence between the intensity and Doppler 
shift variations. The time lag between the two suggests a possible 
7r/2 phase shift, although sometimes the correlation coefficient 
is smaller than 0.3 and the dominant period of intensity is larger 
than that of the Doppler shift. In all of these events, we do not 
see clear signatures of propagating features (inclined stripes) in 
the time-slit diagrams of intensity and Doppler shift. Such an 
observational fact suggests that the phase speed is at least of the 
order of 1200 km s -1 (50" divided by the exposure time 30 s) if 
these oscillations are signatures of propagating features. Such a 
high speed is much larger than the sound speed and comparable 
to the Alfven speed in the corona. 

The event 1 1 (presented in Figure 6(C)) only exhibits oscilla- 
tion in a narrow region (~10") on the slit, which might suggest 
that the loop width is about 10", since the slit was found to make 
a large angle with respect to the loop top (see Figure 8(A)). No 
clear correlation was found for the intensity and Doppler shift 
variations in this event. We have found that the top parts of 
loops associated with oscillating events 12, 13, 14, 15, and 18 
all made a relatively large angle with respect to the EIS slit. The 
coherence scales of these oscillations are typically 10"-30". 
These type II oscillations are usually associated with the quasi- 
periodic appearance of (enhanced) redshift and no significant 
blueshift has been found, which is a puzzle but might be partly 
related to the assumption of zero Doppler shift of the average 
line profile. 

Most of these type II oscillations are found in ARs. However, 
two events (events 10 and 21) are found to be associated with 
coronal BPs. In these two events, we see a relatively large 
intensity fluctuation and a clear 7r/2 phase shift between the 
intensity and Doppler shift oscillations. Analysis results of event 
21 are presented in Figures 6(D)-(F), and the oscillating BP can 
be found in Figure 8(B). The blueshift between y = —20" and 5" 
and redshift between y = —45" and — 20" might be related to 
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Figure 8. Coronal structures associated with two type II oscillation events (events 11 and 21). (A) An XRT image taken at 19:49 on 2007 February 21. The black 
vertical line indicates the location of the EIS slit. The region between the two horizontal bars corresponds to the y- range of the Dopplergram in Figure 6(C). The arrow 
points to the approximate location of oscillation event 1 1. (B) An XRT image taken at 17:53 on 2007 April 19. The region between the two horizontal bars corresponds 
to the y-range of the Dopplergram in Figure 6(D). The arrow points to the approximate location of oscillation event 21. 

(A color version of this figure is available in the online journal.) 


our assumption of zero Doppler shift of the average line profile. 
It could also be caused by a siphon flow along the small BP loop 
system. 

The periods of the type II Doppler shift oscillations are in the 
range of 3-6 minutes, with a mean of 4.6 minutes. The global 
wavelets usually peak sharply and thus the uncertainties of the 
periods are smaller compared to type I oscillations. The aver- 
age oscillation amplitude, as computed from the global wavelet, 
is 1.6 km s -1 . Such a small oscillation amplitude is extremely 
difficult to identify from on-disk imaging observations if they 
correspond to transverse displacement oscillations. With spec- 
troscopic observations, we can clearly see them. 

Unlike type I oscillations, it is difficult to find any significant 
profile asymmetry at the oscillation locations. The almost 
symmetric line profiles seem to exclude the presence of high- 
speed outflows, which are believed to be responsible for the type 
I oscillations. Perhaps a possible mechanism for these type II 
oscillations is the intermittent obscuration by fibrils (De Pontieu 
et al. 2005; de Wijn et al. 2007), or the coronal response of 
fibril motion. Imagining cool chromospheric fibrils going up and 
down quasi-periodically, the overlying hot coronal plasma can 
be forced upward and released downward with the same period. 
Such a scenario seems to be similar to vertical piston motion. In 
this case, the entire coronal plasma is moving at a small speed, 
which may explain the absence of profile asymmetry and line 
broadening. However, it is unknown whether the fibril motion 
can affect the tops of these coronal loops, which are often high 
(e.g., ~50" for the type I oscillation event shown in Figure 1) in 
the corona. 

4.2.1. Mode Identification and Coronal Seismology 

More likely, the type II oscillations are signatures of MHD 
waves. In the following, we discuss possible wave modes 
which might explain our observations and their roles in coronal 
seismology when applicable. 

1. Propagating slow wave. Slow-mode magnetoacoustic 
waves are compressible modes and a correlation (or anti- 


correlation) between intensity and Doppler shift is expected 
in the propagating slow wave. Our observations reveal very 
small intensity fluctuations of all type II oscillation events, 
except events 10 and 21. Moreover, we do not see any clear 
signature of correlation (or anti-correlation) between in- 
tensity and Doppler shift oscillations in any of our type II 
oscillation events. Thus, we can rule out the possibility of 
propagating slow waves. 

The top parts of oscillating loops are more or less parallel 
to the EIS slit in events 17, 19, and 20 (see the example 
in Figure 1). Given a typical coronal sound speed of 
200 km s _1 , the wave would need 180 s (much larger than 
the cadence of 30 s) to propagate through a distance of 
50", so that we should be able to see propagating features 
(inclined stripes) in the time- slit diagrams of intensity and 
Doppler shift. The absence of such propagating features 
again rules out the possibility of propagating slow waves. 

2. Standing slow wave. The seemingly 7t/2 phase shift be- 
tween intensity and Doppler shift in some oscillation events 
suggests the possibility of standing slow waves (Sakurai 
et al. 2002; Kitagawa et al. 2010). For the fundamental 
mode, the wave period is determined by twice the loop 
length divided by sound speed (Wang et al. 2002). Due 
to the strong foreground emission and the mixture with 
emission from surrounding structures, it is often difficult to 
accurately locate the legs of the oscillating loops (so loop re- 
construction techniques are highly desired, e.g., Feng et al. 
2007; Syntelis et al. 2012). Supposing that the oscillating 
loop legs are located at [300", -130"] and [320", -30"] 
in Figure 1 and assuming a semicircular loop, we obtain a 
loop length of 160". Using a sound speed of 200 km s -1 , the 
period was estimated to be 19.22 minutes, which is much 
larger than the observed period (4.91 minutes) of event 20. 
Therefore, the fundamental mode of standing slow waves 
cannot explain type II oscillations in such large loops. A 
possibility that can be considered is a higher harmonic, but 
compared to fundamental mode, harmonic modes are dif- 
ficult to be excited in coronal loops. In addition, for most 
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of the type II oscillations, the correlation between inten- 
sity and Doppler shift is poor (coefficient smaller than 0.3), 
which again does not favor a slow wave interpretation. 

However, standing slow waves may not be excluded for 
oscillations found in smaller loops such as the events of 
February 16 and April 19 (events 10 and 21 in Table 1; 
Figures 6(D)-(F)). In these two events, we also see a 
relatively larger fluctuation of intensity, similar periods of 
intensity and Doppler shift oscillations, and a 7r/2 phase 
shift between intensity and Doppler shift, which favors 
the interpretation of standing slow waves if the intensity 
fluctuations are signatures of density change (e.g., Sakurai 
et al. 2002; Kitagawa et al. 2010). In such small BP loop 
systems, undamped propagating waves from one loop leg 
might reflect when they reach the other loop leg, thus 
producing standing waves. The first oscillation in Erdelyi 
& Taroyan (2008) might be of a similar type. However, as 
we will explain in the following, the intensity fluctuations 
might be related to the fact that different parts of a loop 
are sampled at different times. In this case, the intensity 
fluctuations are not related to density change, but may be 
caused by loop displacement during transverse oscillations. 

3 . Standing kink wave. For the fundamental mode of a standing 
kink wave, the period is determined by twice the loop length 
divided by the kink speed (e.g., Roberts et al. 1984; Wang 
et al. 2002; Aschwanden & Schrijver 2011) defined as 


Ck 


1 + Pel Pi 


Yam 


( 1 ) 


where p, and p e are internal and external densities, respec- 
tively. Vai is the internal Alfven speed, which is related to 
the internal magnetic field (Bp and density in the form of 


Va; = 


B, 

\fVnpi' 


( 2 ) 


In order to explain the period (4.91 minutes) of oscillation 
event 20, the kink speed needs to be about 775 km s _1 . 
This value is in the range (543-2322 km s -1 ) of those 
measured from TRACE transverse oscillations (Ofman & 
Aschwanden 2002). However, we have to bear in mind that 
the value of Ck depends on the magnetic field strength and 
density contrast between the loop and background, both 
of which can vary a lot in different cases. In addition, 
the oscillation could be modulated by some continuous 
driver like p-modes, as suggested by the fact that the 
oscillations are persistent and undamped. Thus, the period 
may differ from the value predicted by twice the loop length 
divided by the kink speed. Also, the above phase speed is 
based on cylindrical or slab geometry approximation of real 
coronal loops. Finally, kink modes can also exhibit some 
intensity disturbances if the angle between the LOS and the 
oscillating structure is not exactly 90° (Cooper et al. 2003; 
Wang et al. 2012b). So the standing kink mode cannot be 
ruled out. The weakly damped Doppler shift oscillation 
event (second event) interpreted as kink mode by Erdelyi 
& Taroyan (2008) may be of a similar type. They did not 
find any flare or prominence eruption associated with the 
oscillation, which is also similar to our type II oscillations. 

If these waves are standing kink waves, then we can in 
principle obtain the internal Alfven speed and magnetic 
field strength according to Equations (1) and (2). Using 


the line ratio method, we have calculated the internal 
and external densities from observations. However, the 
density ratio p e / pi, we obtained is usually in the range 
of 0.2-0. 7, which is much larger than the density ratio 
normally assumed (~0.08). Such results can be explained 
by the fact that LOS summation substantially reduces the 
density contrast (De Moortel & Pascoe 2012). Fortunately, 
the density enters through square root and the dependence 
of the phase speed on the density and the contrast is weak. 
In fact, there is only about 25% difference in Ck between 
the density ratios of 0.08 and 0.7. Here, we simply took 
the typical value of 0.08 (Aschwanden & Schrijver 2011) 
and obtained a value of 570 km s -1 for the internal Alfven 
speed in oscillation event 20. Taking the calculated density 
of log (Ni/ cm -3 ) = 9.14, the magnetic field strength of 
the oscillating loop is estimated to be 11.6 G. Here, the 
internal mass density p t is calculated as 1 .2m p Ni, where 
m p is the proton mass and the constant of 1.2 is due to 
the consideration of helium abundance (e.g., Priest 1982; 
Wang et al. 2009a). Using a density uncertainty of log 
(Aj/cm -3 ) = 0.03, a period uncertainty of 0.67 minutes, 
and assuming a 20% uncertainty of the loop length, the 
uncertainty of the estimated magnetic field (Nakariakov 
& Ofman 2001) would be ~25% (2.9 G). Note that 
Equation (1) is based on the simple configuration of a 
thin and straight cylindrical tube embedded in a uniform 
magnetic atmosphere (Edwin & Roberts 1983; Roberts 
2008), which might be too idealized for our coronal loop 
oscillations. In addition, through numerical simulations, 
De Moortel & Pascoe (2009) have demonstrated that the 
magnetic field derived from coronal seismology might 
differ from the real magnetic field by ~50% in some 
cases. 

4. Propagating Alfvenic wave. Ubiquitous coronal Alfven 
waves have been found by observations of both the Coronal 
Multi-channel Polarimeter (Tomczyk et al. 2007, 2008; 
Tomczyk & McIntosh 2009) and the Atmospheric Imaging 
Assembly (AIA; Lemen et al. 2012) on board the Solar 
Dynamics Observatory ( SDO ; McIntosh et al. 2011). Van 
Doorsselaere et al. (2008) proposed that these waves are 
better called fast-mode kink waves rather than Alfven 
waves. However, another theoretical paper, Goossens et al. 
(2012), presents a different view and clearly states that 
these waves can be considered as surface Alfven waves, 
or Alfvenic waves in the nomenclature of Goossens et al. 
(2009). McIntosh et al. (201 1) also pointed out in the online 
supplementary information that (1) in the highly dynamic 
atmosphere, there are no stable cylindrical waveguides with 
straight magnetic field lines that Van Doorsselaere et al. 
(2008) require for kink waves; (2) the term Alfven wave 
is commonly used by the communities of fusion plasma 
physics and space physics for a largely incompressible 
transverse wave for which the major restoring force is the 
magnetic tension. Here, we follow Goossens et al. (2009) 
and McIntosh et al. (2011) and use the term Alfvenic 
waves to describe more loosely the magnetic oscillations 
we observed with EIS. 

Tomczyk et al. (2007) and Tomczyk & McIntosh (2009) 
found that these transverse waves are most clearly present in 
Doppler shift and that their oscillation period peaks around 
five minutes, which are very similar to most of the type II 
oscillations we found here. Moreover, both these Alfvenic 
waves and our type II oscillations are persistent and no 
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Figure 9. Intensity oscillations associated with Alfvenic waves. (A) Transverse displacement of a loop. X and Y are directions across and along the loop, respectively. 
The region between the two dashed lines (—O'/ 5 ~ 075) indicate the slit location (or targeted pixel, 1" wide). The slit is assumed to be parallel to the loop. The dotted 
line marks the location of the loop center. The color coding is shown in panel (D). (B) Loop displacement relative to X = 2". The period is five minutes. (C) Loop 
velocity. (D) Intensity profile across the loop. The part of the loop between the two dashed lines indicates the slit location at time = 1.25 minutes. (E) Time series 
of the intensity fluctuation and Doppler shift when the angle between the polarization direction and LOS is 60° . The time lag between intensity and Doppler shift 
oscillations and their maximum correlation coefficient are also marked in the panel. 

(An animation and color version of this figure are available in the online journal.) 


significant damping is observed. In addition, the presence 
of the almost vertical strips of Doppler shift features in both 
Figures 5 and 6(B) indicates that the oscillations should 
have a very large propagating speed. Given a phase speed 
of the order of 1200 km s 1 , the wave only needs 30 s 
to propagate through a distance of 50". This duration is 
comparable to the exposure time and observational cadence, 
which can easily explain the vertical strips. All of these 
observational facts suggest that some, if not all, of the 
type II Doppler shift oscillations might be spectroscopic 
signatures of the propagating Alfvenic waves. The periods 
of these oscillations, three to six minutes, suggest that the p- 
modes may greatly modulate the generation or propagation 
of these waves (Tomczyk et al. 2007). The amplitudes of 
these type II Doppler shift oscillations are usually in the 
range of 1-2 km s -1 , which is not that far from the velocity 
amplitude found by McIntosh et al. (2011) in AR using 
AIA data (5 ± 5 km s -1 ). Also, we have to bear in mind 
that the Doppler shift oscillation amplitudes observed by 
EIS are likely to be reduced due to the LOS integration 
and the coarser spatial resolution compared to AIA. For 
a quantification of the effect of spatial resolution on the 
measured Doppler shifts from Alfvenic waves, we refer to 
McIntosh & De Pontieu (2012). 

If these transverse oscillations are propagating Alfvenic 
waves, then we can still calculate the internal magnetic field 
if the Alfven speed is known. Unfortunately, the wave only 
needs about 30 s to propagate through a distance of 50", as 
we mentioned above. Such a short time prevents us from 
deriving the phase speed from the time-distance diagram 
(Tomczyk et al. 2008). 


As we discussed above, various observational facts seem 
to favor the interpretation of Alfvenic waves (kink or Alfven 
waves) for our type II oscillations, although slow-mode standing 
waves could not be excluded for events 10 and 21. 

The Wentzel-Kramers-Brillouin energy flux of these trans- 
verse waves can be estimated as following (e.g., Tomczyk et al. 
2007, 2008; Ofman & Wang 2008; McIntosh et al. 2011): 

F = Pi (Sv) 2 V phase . (3) 

Taking the observed value of the velocity amplitude 
8 v =1.11 km s -1 , the calculated number density of log 
(Aj/cm -3 ) = 9.14, and the derived phase speed of 775 km s _1 , 
we obtain an energy flux of 6.7 x 10 3 erg cm -2 s -1 for event 20, 
which is two to three orders of magnitude lower than the energy 
flux required to balance the radiative and conductive losses of 
the active corona (2 x 10 6 erg cm -2 s -1 ), but as we mentioned 
above, the Doppler shift oscillation amplitudes are likely to be 
greatly reduced due to the coarse spatial resolution and LOS 
integration. This is supported by recent observation showing 
close strands in a loop system oscillating out of phase (Wang 
et al. 2012b). Recent high-resolution observations of SDO/AIA 
clearly reveal an oscillation amplitude of the order of 20 km s -1 
(McIntosh et al. 2011). Such a large amplitude leads to an energy 
flux comparable to that required for heating the quiet corona and 
solar wind. 

4.2.2. Intensity Oscillations Associated with Alfvenic Waves 

Weak intensity oscillations could also occur in the case of 
Alfvenic waves, since periodic loop displacement could lead to 
the scenario that different parts (with different intensity) of a 
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loop are sampled periodically. If the polarization of the swaying 
motion is such that both an LOS component (visible in Doppler 
shifts) and a component in the plane of the sky occur, then 
loops moving into and out of a spatial pixel could lead to slight 
intensity fluctuations that are not a signature of density changes. 

Such a scenario can be easily modeled with a toy model. 
Figure 9 shows how the intensity fluctuation, as well as a 7t/2 
phase shift between intensity and Doppler shift can be produced. 
We assume a loop carrying a simple harmonic Alfvenic wave 
with an amplitude of 3 km s -1 and a period of 5 minutes (panel 
(C)). Such a wave will lead to a swaying motion that causes a 
periodic transverse displacement of the loop by up to ~287 km 
(panel (B)). As a result of this displacement, different parts of 
the loop with different intensities are sampled by the instrument 
slit at different times (panel (A)). Periodic displacement leads to 
periodic sampling of the same part of the loop, thus producing 
periodic variation of the observed intensity. We derived intensity 
profiles across different loops from several raster scans of EIS 
(observed on 2007 February 21, 2007 May 19, and 2007 August 
23) and present in panel (D) a typical normalized intensity 
profile of a loop observed by EIS. Using this intensity profile 
and assuming an angle of 60° between the polarization direction 
( X in panel (A)) and LOS, we obtain an intensity fluctuation of 
~1.5% (panel (E)) if the targeted pixel is ~2" away from the 
loop center. The Doppler shift oscillation has an amplitude of 
~1.5 km s -1 and reveals a clear n /2 phase shift with respect to 
the intensity oscillation. 

Intensity oscillations are usually believed to be signatures of 
compressible waves based on the assumption that the intensity 
change is a reflection of density change (e.g., Kitagawa et al. 
2010). A 7t/2 phase shift between intensity and Doppler shift 
is often believed to be an indicator of standing slow waves 
(e.g., Sakurai et al. 2002; Kitagawa et al. 2010). However, 
from our Figure 9, it is clear that intensity fluctuations might 
also be caused by sampling changing portions of a loop, with 
different densities, due to the periodic transverse motion. In this 
case, a 7r /2 phase shift between intensity and Doppler shift is 
expected if the parts of loop sampled by the instrument slit have 
a monotonically decreasing or increasing intensity profile. Note 
that the n /2 phase shift here is caused by the fact that different 
parts of the loop are sampled at different times and is certainly 
not an intrinsic property of the transverse wave. We emphasize 
that for mode identification when using ElS-like spectroscopic 
observations, the jt/2 phase shift alone is not sufficient and 
one needs to know information such as the associated loop 
structure and phase speed (see Section 4.2.1) for a more definite 
identification. 

The amplitude of intensity oscillation depends on the slope 
of the intensity profile in the parts of the loop sampled by 
the instrument slit, and the angle between the polarization 
direction and LOS. The online movie (ml.mpg) shows how 
these two factors impact the amplitudes of intensity and velocity 
oscillations. Using a fixed angle between polarization and LOS 
(60°), the first part of the movie shows that the intensity 
fluctuation becomes larger when the sampled part of the loop 
has a steeper intensity gradient. Note that the red line marks the 
central location of the sampled parts of the loop. If this sampled 
location is the loop center, then the intensity fluctuation is small 
and the period of the intensity oscillation (if significant) is half 
of the velocity period, since the intensity slightly decreases on 
both sides of the loop center. The second part of the movie 
clearly shows that the intensity amplitude increases as the 
angle between polarization and LOS increases from 0° to 90°. 


Tian et al. 

The velocity amplitude changes in the opposite sense. Except 
for the cases of 0° and 90°, the 7r/2 phase shift is always there. 

As we mentioned above, the intensity oscillation and the tt /2 
phase shift between intensity and Doppler shift variation are 
much more prominent in hotter lines such as Fexiv A274.20 
and Fexv A.284.16 (e.g., Figure 4). This might be explained 
by considering the possible different intensity profiles across 
the loop for different lines. If the average temperature of the 
surrounding coronal plasma is comparable to the formation 
temperature of Fexn A195.12 or Fexm A202.04, then the 
intensity profile should be flatter for these cooler lines and thus 
the possible intensity fluctuation associated with the transverse 
oscillations will be smaller. For hotter lines, there will be a larger 
contrast (steeper intensity profile) between the loop and the 
surrounding (cooler) corona and hence the intensity fluctuation 
will be larger. 

The maximum correlation coefficients are always 0.98 in our 
idealized toy model. In real observations, the sampled parts 
of the loop may not always have a smooth monotonically 
decreasing or increasing intensity profile. This is likely to be the 
case, as a loop is believed to consist of several sub-resolution 
strands. This effect will de- smooth the intensity time series and 
thus reduce the correlation between the intensity and Doppler 
shift. In addition, we assume that the oscillating part of the loop 
is parallel to the instrument slit, which is certainly not always 
the case in observations. If the loop makes a large angle with 
respect to the slit and there is no change of density along the 
loop, then the observed intensity change will be reduced and 
the 7t/2 phase shift might be less clear. In fact, we see a more 
obvious 7r/2 phase shift in events 10, 17, 19, 20, and 21. In all 
of these events, the top part of the oscillating loop is more or 
less parallel to the EIS slit. In other type II oscillation events, 
the oscillating loop makes a large angle with respect to the slit 
and we generally do not see a clear ix /2 phase shift. Another 
major factor that would cause a reduced correlation, a variable 
phase shift and intensity oscillation amplitude is the possible 
change of the sampled loop location throughout the time series. 
This could happen due to either natural evolution of the loop or 
slight wobble/jitter in the EIS pointing. Finally, the instrument 
noise will certainly reduce the correlation coefficient and cause 
a variable phase shift with time. 

5. CONCLUSION 

We have performed a statistical study of Doppler shift 
oscillations using data taken by Hinode/ EIS. We have found 
mainly two types of oscillations: one type (type I) is mainly 
found at loop footpoint regions and the other (type II) is typically 
associated with the upper part of loops. The type I oscillations 
generally show coherent behavior of all line parameters (line 
intensity, Doppler shift, line width, and profile asymmetry), 
apparent blueshift and blueward asymmetry throughout almost 
the entire observational duration. These observational facts seem 
to be consistent with the scenario of episodic high-speed upflows 
(fine-scale recurrent jets). The type II oscillations, with a period 
of three to six minutes, are most clearly seen in the Doppler 
shift and often show no significant variation of the intensity 
and line width. The line profiles do not show any obvious 
asymmetry. These are probably spectroscopic signatures of 
kink/Alfven waves. In addition, we have presented a toy model, 
demonstrating that such transverse waves could also produce 
an observed intensity change and a tt/2 phase shift between 
the intensity and Doppler shift oscillations if the associated 
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loops have monotonic density variation. Thus, the value of 
the phase shift alone is not a reliable diagnostic of the wave 
mode in spectroscopic observations made by instruments such 
as Hinode/ EIS. 
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